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We study the time evolution of two ecosystems in the presence of external noise and 
climatic periodical forcing by a generalized Lotka-Volterra (LV) model. In the first 
ecosystem, composed by two competing species, we find noise induced phenomena such 
as: (i) quasi deterministic oscillations, (ii) stochastic resonance, (iii) noise delayed extinc- 
tion and (iv) spatial patterns. In the second ecosystem, composed by three interacting 
species (one predator and two preys), using a discrete model of the LV equations we 
find that the time evolution of the spatial patterns is strongly dependent on the initial 
conditions of the three species. 

Keywords: Statistical mechanics; population dynamics; noise induced effects; Lotka- 
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1. Introduction 

The study of the effects of noise is a well established subject in several different 
disciplines ranging from physics, to chemistry and to biology i . The essential role 
of the noise in theoretical ecology however has been recently recognized. Some key 
questions in population ecology are related to the understanding of the role that 
noise, climatic forcing and nonlinear interactions between individuals of the same or 
different species play on the dynamics of the ecosystems P]|^]|I||^] |l)ir7llHll??ini1lllllll2j . 
More recently noise induced effects on population dynamics have been investi- 
gated [TTflimiToll 1 611 1 7llT%] . In this paper we study the dynamics of two and three 
interacting species in the presence of multiplicative noise and a periodic driving 
force. Specifically we consider: (a) two competing species and (b) three interact- 
ing species. The multiplicative noise models the interaction between the species 
and the environment, and the driving force mimics the climatic temperature oscil- 
lations. In case (a) the interaction parameter between the species is a stochastic 
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process which obeys an Ito stochastic differential equation. The noise induces spatio- 
temporal patterns, quasi-deterministic oscillations, stochastic resonance (SR) phe- 
nomenon IHIEn] and noise delayed extinction, which is akin of the noise enhanced 
stability [211221 • In case (b) we find noise induced spatial patterns whose time 
evolution is strongly dependent on the initial conditions. Specifically we consider 
two different initial conditions: (i) an uniform initial distribution which gives rise 
to anticorrelated spatio-temporal patterns between the two preys, and (ii) a peaked 
initial distribution which gives rise periodically to a strong correlation between the 
patterns of the three species. 

2. The Model 

2.1. Two competing species 

Time evolution of two competing species is obtained within the formalism of the 
Lotka-Volterra equations |2S1 using the Ito scheme 

dx 

— = fix(a-x- (3 x (t)y)+x£, x (t), (1) 
dy 

— = nv(a-y-0 y (t)x)+yt y (t), (2) 

where £i(t) is Gaussian white noise with zero mean and (^i(t)^j(t')} = a5{t — t')6ij. 
The multiplicative noise models the interaction between the environment and the 
species. The interaction parameters (3 X and (3 y are characterized by a critical value 
corresponding to f3 c = 1. We choose (3 X — (3 y — (3. For (3 < (3 C a coexistence regime 
of the two species is established, while for (3 > (3 C an exclusion regime takes place, 
i.e. in a finite time one of the two species extinguishes. 
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Fig 1. The bistable potential U{0) of the interaction parameter 0{t), centered at /3 = 0.99 
(coexistence region). The values of the parameters are h = 1CT 3 , p = 10~ 2 , n = 0.05. The initial 
value of P{t) is /3(0) = 0.94 (bottom of left well). 

The two regimes correspond to stable states of the deterministic Lotka-Volterra 
model. It is then interesting to investigate the time evolution of the ecosystem for 
(3 varying around the critical value (3 C in the presence of fluctuations, due to the 
significant interaction with the environment. The interaction parameter /3 there- 
fore can be described by a stochastic process which obeys the following differential 
equation 
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df3(t) dU((3) 



dt dp 
where U((3) is a bistable potential (see Fig.l) 



jcos UJ t + 



(3) 



U(0) = h(P - (1 + p)fh A - 2h(fi - (i + P )) 2 h 2 



(4) 



£p(t) is a Gaussian white noise with the usual statistical properties: {£p(t)) = and 
{£p(t)£p(t')) = aS(t — t'). We perform simulations of Eqs. (1) and (JSJ by using a 
bistable potential U(J3) centered at = 0.99 < 1, i.e. in the coexistence region. We 
consider as initial value of the interaction parameter (3(0) = 0.94, which corresponds 
to the left minimum of the potential well (see Fig.l). We fix the intensity of the 
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Fig 2. Time evolution of both population densities x and y, at different levels of the additive noise 
of Eq.(3) and SNR of (x — y) 2 . Species x (solid line) and species y (dotted line). The values of 
the parameters are fi = 1, a = 45, 7 = 10 -5 , u)q/(2tt) = 0.34, = o x = a y = 10 -9 , and the 
initial values of the two species densities are a;(0) = y(0) = 1. (a) For o~a = 10~ 12 no fluctuations 
of the two species densities are observed, (b) For op = 3 • 10 -4 quasi-deterministic anticorrelated 
oscillations appear in the temporal series of the two species densities, (c) For op = 10 — 3 a loss 
of coherence of the temporal series for the two species is observed, (d) SNR corresponding to the 
squared difference of population densities (x — y) 2 . 



multiplicative noise at a constant value a = a x = <r y = 10~ 9 , and we vary the 
intensity of the additive noise. For low noise intensity the two species coexist 
all the time, while by increasing noise one species dominates the other one for a 
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random time and quasi-deterministic oscillations take place. By increasing even 
more the noise intensity the coherent response to environment variations is lost 
(see Figs. 2(a), (b) and (c)). These time variations of competing populations are due 
to noise induced phenomena, namely quasi-deterministic oscillations and stochastic 
resonance. In Fig. 2(d) we report the SNR of the squared difference of population 
densities (x — y) 2 . At a finite noise intensity ap — 5 ■ 10~ 4 the SNR exhibits a 
maximum, which characterizes the SR phenomenon. Now we consider the effect of 
the additive noise erg in Eq. on the extinction time of the species and we neglect 
the multiplicative noise and the periodical forcing. Because of the initial condition, 
the ecosystem is in the coexistence region, then the deterministic extinction time 
of both species is infinity. By introducing noise, exclusion takes place and a finite 
mean extinction time (MET) appears. By increasing the noise intensity we obtain 
a noise delayed extinction: both species survive because of the noise. Then the 
average extinction time grows reaching a saturation value, which corresponds to 
the absence of the barrier of the potential U(f3). The MET as a function of the 
noise intensity erg takes a minimum at a finite noise intensity ap ~ 7.5 • 10~ 4 . We 
perform 1000 realizations of Eqs. J5J and the results are shown in Fig. 3. 
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Fig 3. Mean extinction time of one species as a function of the noise intensity erg of Eq.(3). The 
potential U(fi) and the initial condition for /3(t) are the same of Fig.l. 



3. Spatio- Temporal Patterns 

In this section we study the effects of the noise on the spatial distribution of in- 
teracting species. In particular we analyze the spatial distribution of two different 
ecosystems: (a) two competing species, (b) three interacting species, one predator 
and two preys. To study the spatial effects we consider a discrete time evolution 
model, which is the discrete version of Eqs.(l) and (2) with diffusive terms, namely 
a coupled map lattice [T3"| . 
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3.1. Two competing species 

The time evolution of the spatial distribution for the two species is given by the 
following equations 

<V = ^(l-x^-p^ + ^xl^+D^;-^), (5) 

7 

7 

where x?j and t/"^ are the densities of the species in the site ( at the time step 
n, P n is the interaction parameter at the same time, D is the diffusion constant 
and J2-y indicates the sum over the four nearest neighbors. The X"j and Yfy 
terms are independent Gaussian random variables with zero mean and variance 
unit. The interaction parameter j3 n is a stochastic process obtained by solving 
iteratively Eq.©. The bistable potential is centered at /3 = 1 and /3(0) = 0.95, 
which corresponds to the minimum of the left well. The results of our simulations 
are shown in Fig. 4. 

We observe periodical formation of correlated patterns. After 600 time steps the 
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Fig 4. Spatio-temporal patterns induced by the noise for two competing species. We set /i = 2, 
v = (u;o/27r) = 0.34, a = 10" 5 , a x = a y = 1CT 8 , o = 10" 12 , /9(0) = 0.95, D = 0.05. The values 
of the other parameters are the same used in Fig. 2. The initial spatial distribution is homogeneous 
and equal for the two species, that is a:™' = J/*™ 1 * = 0.5 for all sites 

spatial distributions of the two species appear almost independent. After 3000 time 
steps the anticorrelation between the two species densities increases. A further 
evolution of the system produces again spatial distributions which are independent 
after 4800 time steps and anticorrelated after 5400 time steps. The two competing 
species tend to occupy different spatial places periodically without overlap. 
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3.2. Three interacting species: two preys and one predator 

Finally we analyzed an ecosystem composed by three species: two preys and one 
predator. We use the same coupled map lattice model of the previous section 

= ^[-/3,+7,(^+^)] + ^<^+^E(^-^)' ( 9 ) 

where xfj , • and 2™ • are respectively the densities of preys x, y and of the predator 
z in the site ( at the time step n, 7 and j z are the interaction parameters between 
preys and predator and D is the diffusion coefficient. The interaction parameter [3 
between the two preys is a periodical function whose value, after n time steps, is 
given by 

0(t) = 1 + e + acos(uj t) 7 (10) 

with e = —0.01, a = 0.1 and i/ = (w /27r) = 10~ 3 . We consider two different 
initial conditions: (i) a homogeneous initial distribution and (ii) a peaked initial 
distribution. In the first case we find exactly anticorrelated spatial patterns of 
the two preys, while the spatial patterns of the predator show correlations with 
both the spatial distributions of the preys (see Fig. 5). The preys tend to occupy 




Fig 5. Spatial patterns induced by the noise for three interacting species (two preys and one 
predator) with homogeneous initial distributions. The parameter set is: e = —0.01, fi = 2, 
P z = 0.01, u = (wo/2tt) = 10~ 3 , a = 0.1, cr x = a y = u z = 10~ 8 , D = 0.01,7 = 3 ■ 10~ 2 , 
7z = 2.05 ■ 10 2 . The initial values of the uniform spatial distribution are x 1 ™ 1 = y™ 1 ? 1 = 0.25 and 
z init = g ig for all sitcs 

different positions as in the case of two competing species. In the second case we 
use delta-like initial distributions for the two preys and a homogeneous distribution 
for the predator. After 800 steps we find strongly correlated spatial patterns of the 
preys which almost overlap each other. The maximum of spatial distribution of 
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the predator is just at the boundary of the spatial concentrations of the preys, so 
that the predator surrounds the preys (see Fig. 6). The preys now tend to overlap 
spatially as it occurs in real ecosystems when preys tend to defend themselves 
against the predator attacks. 
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Fig 6. Spatial patterns induced by the noise for three interacting species (two preys and one 
predator) with delta-like initial distributions of the preys and a homogeneous distribution of the 
predator: (a) initial conditions, (b) spatial patterns after 800 time steps. Here we set e = —0.05, 
D = 0.1, o- x = u y = cr z = 10 -3 and the other parameters are the same as in Fig. 5. 



4. Conclusions 

We analyzed the role of the noise on the spatio-temporal behaviors of two ecosys- 
tems composed by two and three interacting species. We use Lotka- Volterra gen- 
eralized equations with multiplicative noise, climatic periodic forcing and random 
interaction parameter. The main conclusions of this work are: (i) in the case of 
two interacting species the appearance of noise induced phenomena such as quasi- 
deterministic oscillations, stochastic resonance, noise delayed extinction and spatial 
patterns; (ii) in the case of three interacting species the formation of dynamical spa- 
tial patterns exhibiting correlations which are strongly dependent on the initial con- 
ditions. Our population dynamical model could be useful to explain spatio-temporal 
experimental data of interacting species strongly interacting with the noisy environ- 
ment 24 ,25 26. A more detailed study of the dynamics of spatial patterns through 
stochastic partial differential equations will be the subject of future investigations. 
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